Gravitational self-force correction to the innermost stable circular orbit of a 

Schwarzschild black hole 
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The innermost stable circular orbit (ISCO) of a test particle around a Schwarzschild black hole 
of mass M has (areal) radius risco ~ QMG/c?. If the particle is endowed with mass M), it 

experiences a gravitational self- force whose conservative piece alters the location of the ISCO. Here 
we calculate the resulting shifts Arisco and Afiisco in the ISCO's radius and frequency, at leading 
order in the mass ratio /x/M. We obtain, in the Lorenz gauge, Ariaco = — 3.269(±0.003)/xG/c^ and 
Af2iaco/fiisco = 0.4870(±0.0006)/i/M. We discuss the implications of our result within the context 
of the extreme-mass-ratio binary inspiral problem. 
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The last few years have seen a breakthrough in the 
development of computational tools for tackhng the 
two-body problem in general relativity Numerical- 
relativistic (NR) codes can now track the complicated 
nonlinear dynamics of a spacetime describing the inspiral 
and merger of two black holes with comparable masses. A 
subset of problems which remains so far beyond the reach 
of full NR treatment concerns the radiative inspiral of a 
small compact object into a black hole of a much larger 
mass. (The study of such extreme-mass-ratio inspirals is 
motivated strongly by the prospects of observing grav- 
itational waves from captures of compact stars by mas- 
sive black holes in galactic nuclei — a particularly exciting 
source for the planned space-based gravitational wave de- 
tector LISA [2].) The large span of Icngthscales and long 
inspiral timescale render this problem intractable using 
current NR technology Q • 

The problem, however, falls naturally within the re- 
mit of black-hole perturbation theory, which utilizes the 
small mass ratio rj as an expansion parameter. In this 
approach, the small object is treated as a structureless 
point particle, which moves on the fixed background of 
the massive hole. At the test-mass limit {rj — + 0) the par- 
ticle traces a geodesic of the black hole geometry. When 
rj is finite (yet small), the particle interacts with its own 
gravitational field, giving rise to a gravitational self force 
(SF). The SF has a dissipative piece which drives the 
radiative inspiral; it also has a conservative component 
which affects the evolution of the orbital phases. 

How to regularize the gravitational self-interaction in 
curved spacetime is now a well understood problem in 
fundamental relativity 3, and a program to calculate 
the gravitational SF in the strong-field exterior of a Kerr 
black hole is well underway. The first actual quantitative 
results were obtained for radial plunge [s^ and circular 
orbits [6] in a Schwarzschild exterior. In the circular or- 
bit case, accurate computations of conservative SF effects 
were performed using two independent methods, success- 
fully tested against each other and against results from 
post-Newtonian (PN) theory [tI, Q . The authors of this 
letter have recently completed the development of a com- 



putational platform for calculations of the SF for generic 
(bound, eccentric) orbits in Schwarzschild. Details of this 
work, which we consider a significant milestone in the SF 
program, are given in a forthcoming paper [9j. 

As a first application of our code, we report here 
the values of the 0{'ri) shifts in the location and fre- 
quency of the innermost stable circular orbit (ISCO) 
of a Schwarzschild black hole due to the conservative 
piece of the gravitational SF. This represents a first con- 
crete calculation of a conservative SF effect with a clear 
physical interpretation and (potentially) observable con- 
sequences. The value of the ISCO frequency, which is 
gauge-invariant (within a class of physically-reasonable 
gauges — see below), has been utilized in the past as a 
convenient reference for comparison and calibration of 
PN schemes [l^. Our result is exact at 0{r]) (within a 
controlled numerical error), and we expect it to provide 
an accurate benchmark against which future PN calcula- 
tions (and potentially also future NR results for binaries 
with small r]) could be tested. 

The radiative decay of a physical circular orbit pro- 
ceeds through an "adiabatic" inspiral epoch, during 
which gravitational radiation reaction slowly removes en- 
ergy and angular momentum from the system, on to 
a brief plunge episode in which the particle drops into 
the black hole along a nearly-geodesic trajectory. The 
transition from inspiral to plunge occurs at radii local- 
ized around the ISCO. The properties of the transition 
regime were studied by Ori and Thorne (OT) [ij] (for 
circular equatorial orbits in Kerr), with later generaliza- 
tions by O'Shaughnessy [11] and Sundararajan [l^. OT 
showed, in particular, that the frequency bandwidth of 
the radiative transition scales as ~ yy^/^: The transition 
becomes more "abrupt" — less gradual — with smaller ry. 
OT (as also [ll,[il) did not tackle the local SF dhectly. 
Instead, they relied on (numerical) computations of the 
flux of energy in the radiated gravitational waves, from 
which the dissipative piece of the SF was inferred indi- 
rectly. OT had to ignore the conservative piece of the SF, 
yet unknown at the time of their analysis. Our work is 
complementary to OT's in that we consider only conser- 
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vative effects and ignore dissipation: it is only then that 
the ISCO becomes precisely locahzable. As we shall con- 
clude below, our result justifies a posteriori the omission 
of conservative SF terms in the OT analysis (at least in 
the special Schwarzschild case). 

In what follows we (i) derive formulae for the conserva- 
tive 0(77) shift in the ISCO radius and frequency in terms 
of the conservative piece of the SF; (ii) describe in brief 
the numerical method used for computing the conserva- 
tive SF (delegating details to a forthcoming paper [§]); 
and (iii) state our results and discuss their implications. 
In the rest of this letter we use standard geometrized 
units, with G = c = 1. We denote by M the mass of 
the background Schwarzschild black hole and by /i the 
mass of the particle (hence /i/M — rj). The particle's 
orbit is described in standard Schwarzschild coordinates 
by i = tp{T), r — rp(r) and ip — (ppir), where r is 
proper time; without loss of generality we take 9p = n/2. 
The metric signature is — I — h-l-, and tensorial indices are 
raised and lowered using the background metric. 

Formula for the ISCO shift: — At the limit /i — > the 
particle traces a timelike geodesic of the Schwarzschild 
geometry, described (using overhats to distinguish the 
geodesic from the perturbed orbit discussed later) by 



tp = E/f{fp), 



= L/fl. 



(1) 



(2) 



Here an overdot denotes d/dr, E and L are, respectively, 
the specific energy and angular momentum (constants 
of the motion), /(r) = 1 — 2A//r, and the radial effec- 
tive potential is VcsirjL"^) — /(I -I- i^/r^). Bound or- 
bits exist only for L > \/V2M, in which case T^ff has 
two extremum points and E"^ — Vdi can have 3 distinct 
roots. In the latter case, the radial motion is bounded as 
r-min <fp{T) < Tniax, where Tmin ( "pcriapsis" ) and r^ax 
("apoapsis") are the second-largest and largest roots, re- 
spectively. For such orbits we define the "eccentricity" e 
and "semi-latus rectum" through 



ro/(l + e), r,, 



ro/(l-e). 



(3) 



The orbits are fully parameterized (up to initial condi- 
tions) by any of the pairs {E, L}, {r^ax, rmin}, or {e, tq}. 
Stable circular orbits occur when E"^ equals the minimum 
of VcS] t hey have e = and radius rmin = '''max = ''o = 
(L2 -I- - 12M2)/(2M). No bound orbits exist for 

L < \/T2M. The circular orbit with L = VT2M is the 
ISCO; it has radius ro = 6M. 

The nature of the ISCO as a separatrix for stable cir- 
cular orbits is understood from the following dynamical 
consideration. Still working at /i — > 0, consider a small-e 
perturbation of a circular geodesic with radius tq, such 
that the resulting slightly-eccentric geodesic has turning 
points rmin and Tmax as in Eq. ^ . For this orbit we write 



where fi(f) is to be determined below. The r component 
of the test-particle's geodesic equation reads 



1 dV.ff 

2 df„ 



= .Foff(fp,L^), 



(5) 



J-'cff being an effective radial force (per unit test mass) 
for the geodesic motion. The linear variation of J-cs 
with respect to e (with fixed ro and f) is de^cS — 
{d!Fcs/dfp)6efp -f [d!Fcs/dL)6eL, where the partial 
derivatives are evaluated at e = 0, and where 5efp = efi. 
5eL is determined by solving E'^ = ^^(rmin, i^) and 
E^ — Vcsifmax, simultaneously for E and i, and then 
using Eqs. ([2]); this gives = Mrl/{rQ - iM - Me^) 
and hence 5eL — 0. At 0(e), therefore, Eq. ([5]) reads 



with 



fi(f) 



,-,2 ^ 



M{ro - 6M) 
r3(ro-3M)' 



(6) 



(7) 



Thus, the 0{e) radial motion is simple-harmonic with 
frequency ujr- The orbit is stable under small-e pertur- 
bations when o)^ > 0, namely for rg > 6M (no circular 
timelike geodesies exist for ro < 3M). The condition 
ujr — identifies the ISCO at rg — 6M. 

We now turn to consider the 0{fi) conservative cor- 
rection to the orbit. The physical SF along the geodesic 
fp(f) is constructed (as prescribed in [^) from the re- 
tarded metric perturbation associated with the particle. 
We denote this SF by fJ-F^^ , and further denote by fJ,F"^^ 
the force derived in just the same way from the advanced 
perturbation [we factor out /i for later convenience, not- 
ing our i^'s(cx fj.) describe self acceleration]. The con- 
servative piece of the SF (per unit /i) is defined through 
F" = (i^.'^t + ^;"dv)/2- It is convenient to take f = 
at a periapsis of rp(f), in which case we have the sym- 
metry relation [TJ] F"^^{t) — €aF^^{~'T) (no summation 
over a), where ea = (—1,1,1,-1) in Schwarzschild coor- 
dinates. This gives us a practical formula for extracting 
the conservative piece from the full (retarded) SF: 



F'^{T) = [FZ,iT) + e^FZ,{-T)] 12. 



(8) 



Notice that, since both F"^ and fp are periodic and even 
in f [and rp(f) is monotonic on < f < 7r/cDr], we may 
express F"^ as a function of fp only (for given ro , e) . We 
write F"^ = F^{fp; ro, e). 

Under the effect of F", the nontrivial components of 
the equation of motion read 



Teffirp, L^)+F\ 



E = -Ft, 



L^F, 



(9) 



(10) 



rp(T) = ro + efi(f) + 0(e2). 



(4) 



where E = tp/ f{rp) and L = ^p/r'^ are (in general) no 
longer constant along the orbit, and hereafter an over- 
dot denotes d/dr. From symmetry, the perturbed orbit 
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remains equatorial. We assume rp(T) remains bounded 
between some Tmin and r„iax, define rg and e as in Eq. 
(I3|), and take r = at a periapsis. Then we may write, 
through O(^), — F''(rp; rp, e). By virtue of Eqs. 
dS]) and (fTO]) we similarly have £' = E{rp;ro,e) and 
i = L(rp; ro, e) through 0(/i). 

Consider first circular orbit solutions of the set (|9ll0p . 
In this case E and L are constant along the orbit. Their 



values are obtained by solving Eqs. |T]) (with 
and Q simultaneously, with rp — r 
This gives 



and Tp = ro. 



rof^iro - 3M)-i [1 - {ro/fo)FS] , (11) 
Mrliro - 3Af)-i [l - (rg/M)F^] , (12) 



where subscripts '0' denote circular-orbit values. 

To identify the new location of the ISCO we again con- 
sider a small-e perturbation of the circular orbit. Writing 
rp(T) as in Eq. we find that r'i(r) again satisfies an 
equation of the form where the frequency is now 



— [^eff(rp,L^(rp))+F^(rp)]^^ 



(13) 



In obtaining this result we have replaced the linear vari- 
ations SeL and S^.F'' with riSr^L and riS^pF^, respec- 
tively, making use of the fact that, through 0(e), L 
and F^ depend on e only through rp(T;ro,e). That 
this is true is shown as follows. Let us formally expand 
L = L(rp;ro,e) = L(o)(rp;ro) + eLW(rp;ro) + 0{e^), 
where the coefficients L*-"-* depend on e only through rp, 
and O(e^) represents terms whose explicit dependence 
on e is at least quadratic. Now note from Eq. ^ that 
formally replacing e — > — e is equivalent to replacing 
Tmin ^ r^iix; hence ^(rmin; ''o, Te) = L(rinax; ro, ±e), 
giving L(i)(rinin;ro) = - L'^^H'Tma^] ro) . Assuming L is 
a continuous function of rp, the last result implies that 
L(i) = 0(e) for aU rp. We thus have L = L(o)(rp;ro) + 
O(e^), with the conclusion that, working at 0(e), L may 
depend on e only through rp . The same argument — with 
the same conclusion — applies to F"^. 

It is now useful to describe the 0{e) orbit more ex- 
plicitly: Assuming periapsis at t = 0, Eq. ^ (over hats 
removed) integrates to give ri = — roCOSw^T, hence 



rp(T) = ro(l — ecosiVrT) + 0{e^). 



(14) 



We can then expand L{rp) through 0{e) in the form 
L = Lo + eriL'(ro) — Lo — eroL'(ro) cosw^t, where a 
prime denotes d/drp and Lq is the circular-orbit value, 
given in Eq. p^ . Since F^p = L, we obtain the form 



Fip — euJrF], sinuJrT + 0{e^ 



(15) 



where F^ — roi'(ro). Similarly, we expand i^''(rp) 
FS + eriF'-'{ro) + Oie^), giving 



F"- = F^ + eFl cosLOrT + 0(6^), 



(16) 



where F[ = —roF'^'{ro)- Replacing L'{ro) F^/ro and 
F'''{ro) -F[/ro in Eq. ([T3D and substituting for Lo 
from Eq. IT^ . we finally obtain 



(17) 



where a = -Srg'^ro - 4Af)/(ro - 3Af), /3 = r^^ and 7 = 
-2rj;"''[Af(ro - 3M)]^/^. This formula describes the O(^) 
conservative shift in the radial frequency off its geodesic 
value. Note it requires knowledge of the SF through 0{e). 

The perturbed ISCO radius risco is now obtained from 
the condition ^^(risco) = 0. Recalling Eqs. d?]) this gives 



risco - 6Af 

{rl/Mmi-ro){aFS 



216M^FS^^ - 108Af2F[i, + V3F'„, 



(18) 



where in the second line we were allowed to replace 
r\sco 6Af since the SF terms are already 0(/i) [so the 
error introduced affects Arisco only at 0{iJ?y\. We have 
denoted F^-^^ = i^o (ro = QM) and similarly for F^^^, F^. 
Note Arisco is gauge dependent, just like the SF itself p^ . 

A more physically-meaningful quantity associated with 
the circular orbit is the azimuthal frequency U, = 
dcpp/dt — ipp/ip. As discussed in [8|, is invariant under 
all 0(/i) gauge transformations whose generators respect 
the helical symmetry of the circular-orbit configuration 
p^ . An expression for the SF-corrected is obtained 
using Eqs. ^ and (fT2| : 



^ _ ro(ro -3A-f) ^^' 
2Af/o °. 



(19) 



where Cl = (Af/rQ)^/^ jg i^^q geodesic (no SF) value. 
Evaluated at ro = rjsco, the SF- induced frequency shift 
Aft = fi — reads [through 0(/i)] 



1 27 

-An,eo/M+-AfF^i3 



(20) 



where ilisco = ^{ro = 6A/) = (6^/2^/)-i Equations 
analo gou s to (fTSl) and (|20p were obtained by Diaz-Rivera 
et al. [I4I in their study of scalar SF effects. 

Numerical calculation of the SF: — Our numerical tech- 
nique for computing the SF is described in detail in 
Rcfs. [6j] (circular orbits) and (eccentric orbits). The 
method is based on a direct integration of the linearized 
(hyperpoblic) set of Einstein equations in the Lorenz 
gauge, with a delta- function source representing the 
stress-energy of the moving particle. The set of 10 cou- 
pled equations is decomposed into tensor- harmonic {I, ni) 
modes, and we solve for each I, m-mode of the (retarded) 
metric perturbation using finite differentiation on a char- 
acteristic mesh in l-|-l-dimensions (time-|-radius). The 
values of the metric perturbation modes and their deriva- 
tives at the particle are then fed into the "mode-sum for- 
mula" [H, [3, 



which yields the physical (retarded) SF 
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through a mode-by mode regularization. The conserva- 
tive piece F" is finally constructed using Eq. ([5]). 

For our ISCO analysis we need the values of the SF 
coefficients Fq , F^ and F^, evaluated at ro = 6M. These 
are extracted by fitting our numerical SF data along the 
orbit using the analytic formulae (fT5|) and (fT6|) . We used 
two independent methods to obtain the necessary data. 
The first is based on applying our eccentric-orbit code [9] 
for a sequence of slightly-eccentric geodesies approaching 
the ISCO. The second method (following the treatment 
of Diaz- Rivera et al. in the scalar-field case [Ol) is based 
on expanding the perturbation equations about a circular 
orbit and solving through 0(e) for a sequence of circular 
geodesies approaching tq = 6Af. In both methods, the 
ISCO limit is quite delicate (as the radial period tends to 
infinity there) and requires a careful treatment — see @ 
for a detailed discussion of our numerical procedure and 
error analysis. Reassuringly, our two methods produce 
mutually consistent numerical data. 

Results and discussion: — We obtained F^-^^ — 
0.0244666 /x/M2, F^-^ = 0.06209 /^/M^ and F\^ = 
— 1.07/x. Equations (fTS]) and (|20|) then give, respectively, 

An^co = -3.269 ^ 0.4870 fi/M. (21) 

^ 'isco 

We estimate the absolute numerical error in these fig- 
ures as ±0.003 /i and ±0.0006 fi/M, respectively. This 
error is controllable and can be reduced in the future 
simply by running our code for longer. Recall the quoted 
values correspond to the Lorenz gauge. Note the ISCO 
frequency increases under the conservative effect of the 
gravitational SF. This is qualitatively similar to the ef- 
fect of the scalar SF acting on a particle endowed with 



scalar charge, which also increases ilisco 

It is interesting to compare the conservative shift 
Ariisco with the frequency bandwidth of the dissipative 
transition across the ISCO. From the analysis of OT [llj 
one estimates the latter (in the Schwarzschild case) as 
ArJdiss/f^isco ^ 4.3874772/5, where, recall, 77 = fj,/M. 
Hence, Afidiss/Afiisco — 977"^/^, giving, for example, 
~ 35830, 9000 and 2261 for mass ratios r] = 10"^, 10"^ 
and 10~^, respectively. Thus, for rj in the astrophysically- 
relevant range, the conservative piece of the SF shifts the 
ISCO frequency by an amount much smaller than the 
entire bandwidth of the transition regime. This confirms 
what one expects based simply on the above rj^^^^ scal- 
ing, and it provides a firm justification for omitting the 
conservative SF in the OT analysis. 

The main practical value of our result, we believe, is 
in that it provides an accurate strong-field benchmark to 
inform the development of approximate methods, such 
as FN. A particularly promising FN treatment of the 
small-?] inspiral problem is based on the Effective One 
Body (FOB) approach , with a recent study ^ con- 
vincingly demonstrating the effectiveness of EOB in de- 
scribing the dissipative part of the strong-field dynamics. 
It would be interesting to test the performance of the 
conservative EOB dynamics against the results obtained 
here. We also envisage that our results could help in- 
form the refinement of approximate "Kludge" models for 
gravitational waveforms from astrophysical inspirals ^2^ . 
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